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The anomalous exponents governing the long wavelength behavior of the flat phase of physical 
crystalline membranes are calculated within a self-consistent screening approximation (SCSA) ap- 
plied to second order expansion in 1/dc (dc is the co-dimension), extending the seminal work of 
Le Doussal and Radzihovsky [Phys. Rev. Lett. 69, 1209 (1992)]. In particular, the bending rigid- 
ity is found to harden algebraically in the long wavelength limit with an exponent rj — 0.789..., 
which is used to extract the elasticity softening exponent r\ u = 0.422..., and the roughness exponent 
£ = 0.605.... The scaling relation rj u — 2 — is proven to hold to all orders in SCSA. Further, apply- 
ing the SCSA to an expansion in 1/dc, is found to be essential, as no solution to the self-consistent 
equations is found in a two bubble level, which is the naive second order expansion. Surprisingly, 
even though the expansion parameter for physical membrane is 1/dc = 1, the SCSA applied to 
second order expansion deviates only slightly from the first order, increasing £ by mere 0.016. This 
supports the high quality of the SCSA for physical crystalline membranes, as well as improves the 
comparison to experiments and numerical simulations of these systems. The prediction of SCSA 
applied to first order expansion for the Poisson ratio is shown to be exact to all orders. 

PACS numbers: 82.45.Mp, 87.16. D-, 61.46.-w, 11.15.Pg. 



I. INTRODUCTION 

Physical membranes are 2-dimensional (2D) surfaces 
embedded in the 3-dimensional space. A subset of these 
are crystalline membranes, known also as tethered or 
polymerized, that have a solidlike structure of 2D lattice, 
usually triangular or hexagonal, with fixed connectivity 

Examples of such systems are plentiful in our world. 
In biological systems, prominent is the intercellular side 
of the plasma membrane of some types of cells, which is 
built of spectrin proteins in a triangular lattice. Different 
examples, in the soft condensed matter field, are mono- 
layers of polymerized phospholipid molecules suspended 
on an air-water interface, using their amphiphilic nature 
[2]. Recently, graphene, a single layer plane of carbon 
atoms, as well as individual crystal planes of other layered 
materials, were isolated experimentally [3JEI], represent- 
ing the ultimate crystalline membrane. Further studies 
have demonstrated the stability of graphene even when it 
is free standing, i.e, without the support of a substrate, 
and tensionless [5]. 

The stability at finite temperature of a flat phase of 
tensionless crystalline membranes may seem like a viola- 
tion of the Mermin- Wagner theorem, which forbids the 
existence of long range order in 2D systems due to di- 
verging thermal vibrations. However, this seeming con- 
tradiction between experiment and theory is resolved by 
introducing out-of-plane fluctuations. These fluctuations 
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induce frustration between the large thermal vibrations 
in 2D and the competing gain in elastic energy, which sta- 
bilizes the globally flat phase even at finite temperatures 

This asymptotically flat phase is perturbed by the 
transverse displacements, whose amplitude diverges with 
the size of the system L, as />, where C is roughness 
exponent. This behavior is a result of the anomalous 
bending energy of the flat phase, that for small wavevec- 
tors q deviates from its constant value and acquires an 
anomalous exponent k_r(q) ~ q~ v , satisfying a scaling 
relation 77 = (4 — D) — 2£, where D is the dimension- 
ality of the surface, i.e. D = 2 for physical mem- 
branes. This picture, as well as other properties of the 
flat phase, is a conclusion of more than two decades of ex- 
tensive theoretical, experimental and numerical research 
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their seminal work, Nelson and Peliti [B] have estimated 
these exponents to be 77 = 1, £ = |, using a one-loop 
expansion of the bending energy, within a self-consistent 
approximation, and assuming a finite renormalization of 
the elastic constants. 

Aronovitz and Lubensky [7] have used an epsilon ex- 
pansion e = 4— D to falsify this assumption, showing that 
the elastic constants obtain an anomalous exponent rj u , 
vanishing at long distances as q Vu . Using Ward identities 
of the rotational group, they achieved a second scaling 
relation r\ u = (4 — D) — 2rj. 

The theoretical evaluation of the anomalous exponents 
which is considered to be the most accurate is within 
a self-consistent screening approximation (SCSA), intro- 
duced in the context of crystalline membranes by Le 
Doussal and Radzihovsky [1]. This approximation is 
found to reproduce known theoretical limits, viz. the 



2 



limit of large embedding dimension d [5] [TU]; the case 
where the surface dimension equals that of the embed- 
ding space; and the leading order e expansion. 

However, when extrapolated to physical membranes, 
the exponents calculated in the different methods vary 
substantially. Extrapolating the large-d expansion leads 
to r) = T) u = £ = 2/3; extrapolating the e expansion 
predicts r\ — 0.96, rj u — 0.08 and £ = 0.52; whereas 
SCSA predicts 77 = —7^ = 0.821..., rj u = 0.358..., 
and £ = 0.590... . A different value, though rather 
close to SCSA, was recently calculated by Kownacki 
and Mouhanna [14 that have applied a non-perturbative 
renormalization group (RG) approach to the problem of 
crystalline membranes, using a simple ansatz for the ac- 
tion, to lowest powers in derivative and field expansion, 
to predict r\ = 0.849. 

Thus, it is important to assess the issue of accuracy 
of the SCSA prediction for physical membranes, in order 
to determine its relevance to this case. However, the 
SCSA is an uncontrolled approximation, as it includes a 
partial summation of infinite amount of diagrams. Thus, 
its uncertainty is unknown. The main goal of the current 
paper is to study the SCSA expansion and its accuracy by 
considering the effect of higher orders of the expansion. 

To date, going beyond the aforementioned approxima- 
tions, and "select" the theoretical approach which gives 
the most accurate results, demanded the use of numeri- 
cal simulations or experimental results. One can extract 
T] u = 0.50(1), C = 0.64(2), and X] = 0.750(5) from the 
Monte-Carlo simulations of Ref. [18], which should be 
compared with 77 = 0.81 and £ = 0.59(2), found numer- 
ically in Ref. [TS] and [17], respectively. Lately [IS], the 
anomalous exponents were calculated for a sheet of car- 
bon atoms (representing graphene without the electronic 
degrees of freedom) , using a realistic carbon-carbon po- 
tential, giving 77 ss 0.85. These results are compatible, 
though not perfectly, with the theoretical approaches of 
SCSA, non perturbative RG and large-d estimates, as 
well as with the existing experimental measurement of 
the static structure factor of the red blood cell cytoskele- 
ton by small-angle x-ray and light scattering, yielding a 
roughness exponent £ = 0.65(10) [T6"| . 

The paper is built as follows. After a description of 
the model and its long-wavelength properties, it is shown 
that the scaling relation r\ u = 2 — 2r\ is exact to all or- 
ders of SCSA in D=2. Then, the second order SCSA 
is developed. Explicitly, it is demonstrated that a two 
loop expansion, which is a naive second order, has no 
solutions for the SCSA equations. A solution, however, 
appears when using a second order in l/(d — D), lead- 
ing to rj = 0.789... for physical crystalline membranes. 
The paper ends with a discussion of the meaning of the 
results. 



II. FLAT PHASE OF A CRYSTALLINE 
MEMBRANE 

In general dimensionality, i.e., a D-dimensional mem- 
brane embedded in a e?-dimensional world, a useful 
parameterization of a membrane is the Monge rep- 
resentation. Displacements inside the membrane are 
parametrized using a /^-dimensional phonon field u, and 
the out-of-plane field by a dc — d — D dimensional field 

h. As a result, a particle located on the unperturbed 
(flat) membrane in coordinate x, is displaced due to the 
perturbations to a location r = (x + u,h) (in d dimen- 
sions) . 

One assumes an asymptotically flat geometry with 
small out-of-plane perturbations, such that u and h are 
functions of x. To leading order in field gradients, the 
resulting free energy is a sum of the membrane's elastic 
and bending energies: 

F[u, h \ = \f d °x [k(V 2 /i) 2 + 2/*«««*j + Xu l] • C 1 ) 

Here, k is the bending energy, A is the first lame con- 
stant, /i is the shear modulus, and is the strain 
tensor, defined by Uij = \ (d{Uj + djUi + dih ■ djh), for 

i, j = 1,...,D. The in-plane phonon fields appear only 
quadratically, thus can be integrated out of the free en- 
ergy. The resulting effective energy depends only on h, 
and in Fourier space receives the form [T]: 

Feff[h] = l J ^{ K q*\h q f+ (2) 

f d D k f d D k' R^){k,k',g) u \ 
+ J (2ir)° J (2ir)D 4(d - D) k ' n ?-Z n k' n -t-k> j • 

The integrated out elasticity is hidden in the effective 
four point coupling, 

Rt D \k,k',q) = 2 M (fcP T (« ? -)fc') 2 + ^ Y (feP T ®fc)(fc'P T ('?1fc'), 

(3) 

where the transverse projection operator is defined as 
Pj(g) = \5ij — ^r^j- To this end, two main results 
should be emphasized: 

• The dependence of the effective interaction on the 
dimension of the embedding space is trivial, and 
goes like 1/dc- As a result, every interaction of 
fom-h fields will contribute a factor of 1/dc to the 
calculated observable. 

• The h field propagator has the property 
(h(q)h(-q)) = I dcXdc k B TG(q), where ld G xdo 
is the identity matrix in dc dimensions [ks is 
Boltzmann constant and T is the temperature). 
Thus, in the diagrammatic formulation of the 
theory, any closed loop of h fields propagator will 
contribute a factor of dc to the diagram. 
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III. SELF CONSISTENT SCREENING 
APPROXIMATION 



l3K u 




g — k 



FIG. 1: Feynman rules for the physical theory. 



In view of these points and Eqs. ([2] -[3]), two perturba- 
tive procedures to solve the system are apparent. The 
first is to organize the diagrams in powers of l/dc- Ev- 
idently, a diagram with Lh loops and Nr interactions 
will contribute to the nth order in this expansion iff 
n = Nr — Lh- The other perturbative expansion is in 
the number of loops, which basically implies that the ef- 
fective four point elastic interaction is small compared to 
the bending energy. However, this criterion is valid only 
for length scales smaller than J n 2 /ksTY , while we are 
interested in the long wavelength behavior of the system 
(Y is some combination of the elastic constants). 

Let us now concentrate in the physical case of D = 2 
and d — 3. In this case the interaction is completely 
separable: R^(k,k' ,q) = K [qxk] 2 [qx k'] 2 , where K — 
2^fX 1S ^ ne ^ Young modulus, and q = q/q. The 
separability allows writing the Feynman rules of Fig. [T] 
for the h field and for the screening of the interaction. 
The scale evolution of these fields is governed by the two 
point diagrams. Thus, we define: 



SCSA is an extension of a consistent perturbative ex- 
pansion, in which one replaces every propagator with 
the dressed propagator, and every interaction with the 
screened interaction, however cutting the series of dia- 
grams taken into account in the calculation of £ and *Sf 
according to the order of the original expansion. Then, 
one solves Eqs. (|5] - [6j) self-consistently. For example, in 
the seminal work of Radzihovsky and Le Doussal pQ , they 
extended a leading order expansion in l/dc of £ and 'J: 



*(<?) 



d 2 k 
d 2 k 



z(qk) 



[k x q] 4 



y(q\q-k\)\q- 

[k x qf 



(7) 



{^) 2 y(qk)y(q\q~k\)\q-k\^ 



Note that we use an integration variable scaled by q, mak- 
ing it easier to analyze the momentum scaling of the the- 
ory, and its long wavelength behavior. They then solved 
these equations self-consistently in the long-wavelength 
limit q — > 0, assuming an algebraic behavior for the cou- 
plings, and defining, 



h(a,b) 



d 2 k fcn [k x g] 4 
(2tt) 2 \q-k\4-b 

3 r(i-gffir(i + f)r(i + f) 
f)r(2-|)' 



(8) 



167rr(2- 



«±6)r(2- 



Onc notices that the constant in Eqs. ([5] -[6]) can be ne- 
glected when q — > 0, so the SCSA equations in the long- 
wavelength limit are: 



Kn(q) 



*(?) = 



K R (q) 
K 



(4) 



with K R (q) = k B T{(h(q)h(-q))q 4 )-\ and K R (q) the 
screened interaction. This scale evolution can be writ- 
ten in terms of Dyson equations [20 : 



y(q) 



(5) 
(6) 



Here, q 



K k B T 



E(g) is the sum of all 1PI two- 



point diagrams, and ^(q) is the sum of all 1PI four- 
point diagrams. The latter have a symmetry factor of 
| which is factored out of ^(q)- In the long- wavelength 
limit one assumes a critical behavior for the couplings, 
V(q) = yo(q/qo)~ n and z(q) = z (q/q )^, which imply 
the long wavelength behavior to the 1PI amplitudes. 



yo 



2-V-Vv 



zoVo 1 h(r)u,il) 



,qo 



-Vu i /_ \2-2r, 



In order for these equations to be valid for all q (in the rel- 
evant regime of q <C qo), the scaling relation r\ u = 2 — 2n 
should be fulfilled. Dividing the equations one reaches an 
equation for rj, which can be solved analytically to give 
Radzihovsky and Le Doussal have solved for 



r\ at general D and d, and showed that this approach re- 
covers other approximations, viz large d limit, expansion 
in small e = 4 — D, and d = D. The latter is unique to 
the SCSA for crystalline membranes, as it is not the case 
for the SCSA for O(n) model [H]. 

Albeit these advantages, the SCSA is an uncontrolled 
approximation, as it includes a partial summation of di- 
agrams mixing different orders. Thus, it is useful to go 
beyond the one-loop approximation and check deviations 
from it. 
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IV. HIGHER ORDERS CONTRIBUTION TO 
SCSA 

First, let us comment regarding the effect of including 
higher orders on one of the aforementioned results of the 
SCSA - the prediction of the scaling relation r] u — 2 — 2r\. 
It can be easily demonstrated that this result survives all 
orders. This can be concluded by analyzing the topol- 
ogy of the diagrams. Let us assume a diagram of specific 
order. Creating a higher order diagram from this dia- 
gram can be achieved in two ways (within SCSA): (i) in- 
serting an interaction connecting two h field propagators 
(solid lines); (ii) explicitly plotting one of the diagrams 
resummed in lower SCSA orders, and repeating point (i) 
for this new diagram. 

Clearly, only point (i) has to be analyzed, as lower 
SCSA orders are assumed to fulfill the scaling relation. 
Let us assume that the lines we cut carry momenta qk 
and qk' . We cut these by an interaction with momen- 
tum qk, over which we integrate (so that the integration 

measure is q 2 d 2 k, and the integration variable is dimen- 
sionless) . Following the Feynman rules, the new diagram 
will differ from its ancestor by the additional factor: 



d 2 k 



{-f}K R {qk))[k X k] 2 [k X k 



'12 



( 2?r ) 2 0K R {q\k + k\)\k + fc| 4 /3/c H (g|fc - k'\)\k - k'\ A 

(10) 

Taking the long wavelength limit, and moving to the di- 
mensionlcss elasticity and rigidity: 



-OoZ/o 2 ) 



d 2 k 



-2+r)„+2r; 



[jfc X k] 2 [k X k' 



(11) 



(2tt) 2 /c' 4 ->Hfc + £| 4 - ? '|fc- fc'| 4 - ? '' 



In order for this addition not to affect the q behavior, one 
has to demand that the power will vanish, thus reproduc- 
ing the scaling relation. This can be easily checked to be 
valid at general dimensionality 

Additional conclusion is that the general form of the 
SCSA equations at the long wavelength limit are: 



1 = z o Vq 2 cr(z Q y Q 2 ;tj), 
1 = z yQ 2 tp(z yQ 2 ;i]). 



(12) 



(?) 



2-tj 



Uo 



and ip 



the momen- 



Where a = 

If^-J Vo^il)' are independent of 
turn g, and depend on the universal amplitude z^y^ 2 
and critical exponent rj. Moreover, a (ip) is a polynomial 
in zay^ 2 , whose coefficients have alternating signs and 
are functions only of r\. The power of the polynomial 
equals the number of internal interaction lines minus 
one (internal interaction lines) taken into account in the 
calculation of a (-0). 



For completeness, though it is a deviation from the 
physical case D = 2, we discuss one of the most impor- 
tant results of SCSA, predicting that crystalline mem- 
branes are auxetic, i.e., have a negative poisson ra- 
tio. In particular SCSA predicts, for all D ^ {1,2}, 
lim„^a /l = — rn~r • Here, we show that this result 
can be extended to the Dyson equation level, i.e., to all 
orders of SCSA. The proof follows the lines of Ref. pQ. 
We start by rewriting the general dimensionality interac- 
tion of Eq.jji]) as R^ D \k, k', q) = k i kjk' l k' m p(q) i jj m , with 
p(q) = 2fiM(q)+2bN(q), and b = p{2p + D\)/(2(i + A). 
M and iV are defined as (assuming D ^ 1): 



D 

1 , 



-p£(£)pL(£> 



Pjl(q)Pj m (q) +Pl n (q)Pj l (q)) - N ijM 



We note that M and N are orthogonal for all D ^ 2 
(since M(D = 2) = 0). We proceed by writing a 
Dyson equation for the interaction, generalizing Eq. (J6j> , 

as 'pr = *p — *p ■ II • *p*R, where II is the four point IP I 
amplitude. Due to the projection operators II = STr(q), 
where S is the totally symmetric tensor Sijj m = 8ij t i m + 

fiu,jm + $im,ji- Using the orthogonality of M and ./V for 
D =/= 2, one can now write a Dyson equation for each of 
the elastic coefficients /i and b: 



b R (q) 



l + (D + l)fm(q) 



1 +2^7T( 9 )' 



Taking the long wavelength limit, n(q) 
the universal result lim„ 



->u, 



one gets 



for D ^ {1,2}. 



Mg) - ?_ 

'9^0 n(q) - D+l 

As this result does not depend on the order of the expan- 
sion, I conjecture that this is the exact result. 

Alas, this proof cannot be extended to D — 2 due to 
the separability of the interaction in this dimensionality. 
Monte-Carlo simulations have shown that this result ex- 
tends to D = 2, giving v = —0.32(4), compared to the 
v = —1/3 SCSA extrapolation 22], indicating that an 
analytic continuation of this prediction to D = 2 is ex- 
act. A different result was calculated in Ref. [T7]. There, 
a molecular dynamics approach led to a long-wavelength 
Poisson ration of v — —0.16(1). Hence, more numerical 
investigations are needed to settle this difference. 



V. SCSA FOR PHYSICAL MEMBRANES 
APPLIED TO SECOND ORDER 

As above mentioned extending the calculation to sec- 
ond order is not a fully defined procedure. In this section 
we will study this issue, and check the exponents implied 
by each procedure. In Fig. [2] the SCSA expansion is pre- 
sented in a diagrammatic form. If the second order is just 
a two-bubble expansion, one has to cut these sums and 
neglect the last diagram in each sum. In order to use a 



5 



— n— <T)—n\. 

O CD- o> 



■= — +~ • 



(<0 



FIG. 2: Diagrammatic form of the self consistent equations 
for the (a) h field propagator and (b) interaction. A solid line 
corresponds to the h field propagator, whereas a dashed line 
corresponds to the interaction. Thick lines are dressed. The 
equations shown here correspond to the SCSA applied to sec- 
ond order expansion in 1/dc- In the first order SCSA [I] one 
takes into account these equations only up to the one-loop di- 
agram (the first two diagrams in each equation). A naive two 
bubble expansion takes into account the first three diagrams 
in each equation, (c) Second order interaction vertices. 



1/dc expansion, one has to include also those diagrams, 
which include more than two loops, however satisfying 
the condition 2 = Nr — L^. A priori only the 1/dc ex- 
pansion is guaranteed to lead to a correction to a lower 
order theory, as only the 1/dc based SCSA reduces to 
a controlled method. As will be shown, this fact reveals 
itself in the case of the flat phase of physical membranes, 
as no solution is found to the two-bubble expansion. 

In an algebraic form, the SCSA equations to order 



l/d c are (x = z y ): 



ai(r])x 
bi(r))x - 



- a 2 (?7)a; 2 
b 2 {ri)x 2 - 



- a 3 (?7)x 3 , 
b 3 (v)x 3 . 



(13) 



In a two loop expansion one cuts these equations at two 
coefficients, i.e., taking a 3 = 63 = 0. To put the coeffi- 
cients in an algebraic form, it is convenient to define the 
following integral, representing the second order correc- 
tion to the interaction vertex, i.e. an internal diagram in 
which three legs are connected via a loop (see Fig. [2jc)): 



this definition the coefficients are: 
ai = Ii(Vu,v), 



I{a, b, c; k) 



d 2 k' \k' x fc| 2 |fc' x g| 2 |(fc' - k) x (k' - g)\< 
(2tt) 2 k> 4 - a \k' - k\ 4 - b \k' - q\ 4 - c 



(14) 

where for abbreviation we used the notation k = (k,9), 
and q = (1,0). Evidently, in our case only the integrals 
with at least two of (a, 6, c) equal 77 contribute. Using 



n 2 



d 2 k [k x q} 2 
(2tt) 2 k 4 -v\q - k\ 4 -*h 
! f d 2 k [k x q} 2 



a>3 



b, = 



2 J (2tt) 2 k 4 -n\q - k\ 4 -n 

d 2 k I(r), T), 77; k) I(t) — 2, r], 77; k) 
(2tt) 2 k 4 ~^\q-k\4-v u 



I(r) u ,r),r);k), (15) 



1 f d 2 k 



2 J (2^) 2 fc4-„ u |^_fc|4-„, 



-I 2 (r),r],r}; k). 



VI. RESULTS AND DISCUSSION 

To accomplish the calculation of coefficients 01, . . . , 03 
and b±, . . . , 63, I use a two dimensional Gauss quadrature. 
The calculation starts with a tabulation of the function 
I(a, b, c; k), which takes 0(N A ) operations, with being 
the number of integration points in the range (0, 1). All 
other integrations are 0(N 2 ), thus the tabulation deter- 
mines the length and scaling of the calculation with N. 
The coefficients are calculated for each 77, and are used to 
solve each of Eqs. (13 1 for zqi/q 2 . If the latter are equal, 
then 77 is a solution to the SCSA equation. In order to 
reach a result convergent to 10~ 3 it is sufficient to use 
N = 200 [23]. 

The calculation shows that when considering a two- 
bubble expansion, one cannot find a solution to the SCSA 
equations, i.e., Eqs. (13 1 with 03 = 63 = 0. However, a 
calculation to second order in 1/dc leads to 

77 = 0.78922(5), 

as well as zoUq 2 = 12.763(5). The latter can be put as a 
universal relation between the bending rigidity and the 
Young modulus at large wavelengths: 



lim 

q ^a q 



1 k B TK R {q) 



= 3.573(1). 



When compared to first order SCSA [T], these results 
give a measure of the accuracy of the SCSA. In first order 



SCSA, 77 



0.82085... and 2 % = 11.2276. 



1+V15 

Evidently, these are close in values to the second order 
SCSA. In fact, the critical exponent 77 changes only by 
about 0.03. Together with the fact that the SCSA for 
crystalline membranes coincides with other theoretical 
estimates in the limits dc — >• 00, small e = 4 — D, and 
dc — 0, this is a signature for the relevance of this ap- 
proximation for the physical case d — 3, and D = 2. 

The origins of this success are unclear, and even en- 
hance when considering how essential the description of 
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SCSA as a 1/dc expansion, i.e., the expansion parameter 
is 1, as was shown here. One explanation might be hidden 
in the fact that contrary to SCSA of other field theories, 
the most famous is the 0(n) model [2T], the correction 
to the interaction vertex is finite, and does not demand 
any regularization. Additional sources might be hidden 
in the symmetries of the D — 2 problem [24J . 

The current calculation predicts the other two scal- 
ing exponents as well, r\ u — 2 — 2r\ = 0.4216(1), and 
C = 0.60539(3). The latter is in very good agreement 
with the experimental value £ = 0.65(10) 16J, as well as 
with the molecular dynamics value ( — 0.59(2), and goes 
in the correct direction to the Monte-Carlo simulation 
result £ = 0.64(2). The calculated rj agree, on average, 
with the different numerical simulations, whose results 
are in the range r) — 0.75 — 0.85 [T5J [T5J QI5]. The current 
calculation is also compatible with the non-perturbative 
RG result rj = 0.849 [14], considering the fact that the 
latter included only low powers in the derivative and field 
expansion for the action. Stability is found in this ap- 
proach as well, as the authors find the result stable to 
changes in the cutoff function. Extending these two non- 
perturbative approaches to pinpoint the theoretical pre- 
dictions for the properties of physical crystalline mem- 
branes, is highly called for. A rigorous way to estimate 



the quality of SCSA is by computing the flat phase fixed 
point to order e 2 , and comparing it to the SCSA pre- 
diction of the second order e expansion. The difference 
between the predictions of these theories, if found, will 
indicate the level of precision of SCSA [25] . 

It would be interesting to check the current work us- 
ing experiments and numerical simulations for various 
interparticle potentials, intended to probe the universal- 
ity of the membrane theory, not only for the critical ex- 
ponents, but also for the universal relation predicted to 
exist between the bending rigidity and the Young mod- 
ulus at large distances. Of particular interest is to verify 
whether the exact result for the Poisson ratio at long 
wavelength that was given here for all membrane dimen- 
sionality, excluding D = 1, 2, can be extended to these 
physical cases. 
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